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1 Abstract 

<N . 

Thanks to post-flight analyses of several artificial satellites carried 
out during last years, the meteoroids hazard for space navigation and 
OO ' in-orbit satellites permanence is now clear. Even if catastrophic im- 

pact is a rare event, high meteoroids fluxes can erode and weaken the 
satellite or space station main structures. However, the main danger 
• | seems to be the impact-generated plasma, which can produce electro- 

^ magnetic interferences, disturbing the on-board electronics. 

Jl ' 

PACS-96.50.Kr Meteors and meteoroids. 
PACS-95.40.+S Artificial Earth satellites. 
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1 Introduction 



Till now, it was a common idea that space debris were the main threat 
for spacecrafts. This hypotesis finds its pillar in a paper by Laurance and 
Brownlee (1986). The work was the result of the Solar Max Satellite post- 
flight analysis, from which was evident that space debris flux was several 
hundred times higher than natural meteoroids flux. 

However, during last years emerge some facts that claim for a different 
reality. During 1993 Perseids maximum, Beech and Brown (1993) launched 
an alarm, because they noted that impact probabilities with Hubble Space 
Telescope sized objects, were low, but not negligible. Really, on the night of 
August 12, 1993, astronauts onboard the Mir Space Station reported audible 
meteoroids impacts and, then, it was verified that Mir experienced about 
2000 hits during 24 hours and solar panels were hardly damaged (Beech 
et al, 1995). In the same night, the ESA (European Space Agency) lost the 
control of Olympus telecommunication satellite. The following investigations 
made clear that the failure was probably caused by an impact with a Perseid 
meteoroid (Caswell et al, 1995). 

The satellites LDEF (Long Duration Exposure Facility) and EURECA 
(EUropean REtrievable CArrier) take the lion's share in these studies. Post- 
flight analyses on these satellites drastically rescaled the theories of Laurance 
and Brownlee on the space debris danger. McDonnell et al. (1997b) showed 
that the error in the paper of Laurance and Brownlee was due to the use of 
a non correct formula for transformation from crater dimension to particle 
mass. Then, this scientists group elaborated a new formula, by using lab- 
oratory simulations (Gardner et al., 1997). New fluxes estimates, based on 
data collected on EURECA and LDEF, showed that, at micrometer dimen- 
sions and 500 km altitude, debris population is not dominant as previously 
thought (McDonnell et a/., 1997b). Really, above 30 /iin ballistic limit me- 
teoroids dominates, while, for thickness between 4 and 5 /im, 18% only of 
impacts are due to interplanetary matter. Moreover, the debris flux does 
not change appreciably in the 30 /im size regime over the period 1980-1994, 
owing to atmospheric drag. 

It is then very important to know the distribution and dynamics of inter- 
planetary matter, in order to properly plan satellite orbits and space probes 
courses. 
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2 Distribution and dynamics of interplane- 
tary matter 

In interplanetary space, there are great matter quantities, coming from dis- 
ruption of asteroids and comets. The mass range is very wideQ and is observ- 
able by using several techniques: visual, radar, in situ and others. Crossed 
references among data obtained by these techniques are full of difficulties, 
but not impossible (Ceplecha, 1992; Foschini, 1997). 

Interplanetary matter is subjected to gravitational fields of Sun and plan- 
ets. In a certain sense, it is possible to say that planets are the "road-sweeper" 
of the Solar System, because they collect all these particles that are in their 
neighbourhood for a sufficient time, making some gaps in meteoroids spatial 
distribution (Opik, 1951; see also the figure at p. 240 in Lindblad, 1987). On 
the other hand, the solar radiation pressure, due to the momentum carried 
by solar photons, pushes meteoroids toward outer space. In a heliocentric 
reference frame, let r be the radial unit vector and let 6 be the unit vector 
normal to f in the orbit plane; the meteoroid speed is then: 

v = rf+ r66 (1) 

and the radiation force can be written as: 



mv ~ Q pr ( — ) 
c 



c c 



(2) 



where SA is the total amount of energy intercepted, per second, from a radi- 
ation beam of integrated flux density S by a stationary, perfectly absorbing 
meteoroid of cross section A and Q pr is the radiation pressure coefficient, pro- 
portional to the total momentum withdrawn from the beam (Burns, 1987). 
Usually, the velocity-dependent part of (2) is called Poynting-Robertson ef- 
fect, while the radial term is simply called radiation pressure, even if, as 
Burns (1987) wrote, there are also other accepted customs. 

There is another non-gravitational force, the so-called Yarkovsky effect. 
This effect is connected with the name of a polish engineer who first described 
it in a pamphlet published in russian around 1900 (see Opik, 1951). It is 

1 Generally speaking, a 'meteoroid' is intended to be a cometary or asteroidal body 
with a mass range between 10~ 9 and 10 7 kg, even if there is a vivid debate around this 
definition. 
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a force generated by asymmetries in the reradiated thermal energy from 
a rotating body exposed to Sun, because the evening hemisphere is slightly 
warmer than the morning one. The warmer hemisphere radiated more energy, 
and hence momentum, than the other hemisphere, producing a net force, 
which depends on rotation frequency, thermal properties and dimensions of 
the body, and from Sun distance (Burns, 1987). In a similar way, it is 
possible to speak about a "seasonal" effect, where the temperature difference 
between summer and winter produces a net force, which depends on polar 
axis orientation (Rubincam, 1995). If collisions change the rotational state of 
cosmic bodies, the amplitude and direction of the force due to the Yarkovsky 
effect change in time in a stochastic way, producing a random walk of the 
major semiaxis of the body orbit. 

The ratio among these forces can vary according to dimensions and masses 
of interplanetary bodies, distributing the matter in the Solar System. Only 
those particles, that are in the neighbourhood of a planet for a little time 
thanks to Poynting-Robertson effect, can survive to the gravitational attrac- 
tion ("road-sweeper" effect). The Poynting-Robertson effect become funda- 
mental for bodies in the range from 0.1 /im to some centimetre (Burns, 1987). 
The Yarkovsky effect is dominant for larger bodies (0.1-100 m). In the aster- 
oidal population, this effect produces an orbital major semiaxis drift, driving 
the delivery of meteoroids toward Earth (Farinella et al, 1998). 

Taking into account observations, made with several techniques, and the 
dynamics of the interplanetary matter, it is possible to elaborate a model, in 
order to know meteoroids fluxes and concetrations. Now, the starting point is 
the model elaborated by Griin et al. (1985) and regarding the interplanetary 
matter dynamics at 1 AU (Astronomical Unit) from the Sun. Throughout the 
model, the following characteristics of meteoroids population are considered: 
(1) the mean mass density is 2500 kg/m 3 ; (2) the relative speed between dif- 
ferent meteoroids, as well as the impact speed on the Moon, is 20 km/s; (3) 
the flux onto Earth, as well as mutual impact flux, is isotropic. Divine (1993) 
extended this model, taking also into account data from space probes (Pio- 
neer 10 and 11, Ulysses, Galileo and Helios-1), and found five populations of 
meteoroids, each named for a distinctive characteristic of their orbital distri- 
butions and of their mass. Taylor and McBride (1997) used data recorded 
by Harvard Radio Meteor Project in order to extend the model by Griin et 
al. (1985), taking into account the anisotropies in the meteoroids environ- 
ment. This radiant-resolved meteoroid model was precedeed by a reanalysis 
and correction of speed distribution obtained from radars, which can affect 
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large meteoroids distribution in Divine's model (Taylor, 1995). However, 
Griin et al. (1997) recently developed a new model from the Divine's one, 
including the effect of radiation pressure on the meteoroids speed and consid- 
ering impact directions and speeds. The new model lead to four meteoroids 
populations on elliptical orbits and one moving on hyperbolic orbits. 

At last, it should be pointed out that other models are now available 
(Gor'kavyi et al, 1997; Wasbauer et a/., 1997). 

3 Impact probabilities 

The higher risk regions are those related to meteoroid streams, that when 
encounter the Earth give rise to meteor showers. In these streams, the me- 
teoroids number per volume unit can reach very high values, with geocentric 
speeds up to 71 km/s (Kresak, 1993; Jenniskens, 1995). By the time, it 
is known the last meteor storm generated during 1966 by Leonids, when a 
ZHR {Zenithal Hourly Rate) of 150,000 was recorded^. The next perihelion 
passage of Leonid parent body, the comet P/Tempel-Tuttle, happened on 
February 1998. Then a new meteor storm is expected. If it will be so, there 
will be a considerable hazard for artificial satellites (Beech and Brown, 1994; 
Foschini and Cevolani, 1997). 

A rough evaluation, of the impact probability of a meteoroid with an 
artificial satellite, can be made assuming an uniform and constant flux inter- 
cepting an area A [m 2 ] for a time t [s] : 

/ = nFVAt ■ 10~ 13 (3) 

where V [km/s] is the geocentric mean meteoroids speed and n [km -3 ] is the 
spatial number density of meteoroids contained, during normal conditions, 
into an equivalent 1000 km sized cube. During storms or outbursts, it is 
necessary to introduce an enhancement factor F, that for Leonids ranges 
from 300 to 10000 (Beech and Brown, 1994). In order to calculate the spatial 
number density, it is necessary to know the flux density of meteoroids $ 

2 Jenniskens (1995) rose some doubts on this ZHR evaluation, that do not find a 
consideration in radar data. According to Jenniskens, the correct value should be 
ZHR=15,000±3,000. 
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[km 2 -h x ] in a specific mass range, because n is denned as (Koschack and 
Rendtel, 1990a, b): 

$ 

n = (4) 

3600F V ; 

The flux evaluation is very difficult, because the Earth encounter a part 
of meteoroids only. Moreover, it must be take into account the origin of me- 
teoroids: when the parent body is a comet, a part of the nucleus surface only 
is active^ and, then, particles are expelled along several preferred directions, 
arranging themselves on slightly different orbits, even if centered along the 
orbit of the parent body. Furthermore, considering the forces described in 
the previously section, meteoroids will move along particular configurations, 
composed by several filamentsQ. The effect on radar and visual observations 
of these structures is that there are one or more peaks during shower activ- 
ity. Generally speaking, flux evaluation are made by using visual observations 
and taking into account several corrective factors, in order to offset for limit- 
ing elements (Koschack and Rendtel, 1990a, b; Jenniskens, 1994; Brown and 
Rendtel, 1996). However, an important limit is due to maximum magnitude 
(+6.5), corresponding to a meteoroid of about 10~ 5 kg (Jenniskens, 1994). 
But with radars, it is possible to reach magnitude +16, that should corre- 
spond to a mass of about 10~ 9 kg (Zhou and Kelley, 1997). Other advantages 
are a wide collection area (for forward-scatter radars it is possible to reach 
some thousands of square kilometre) and the absence of weather, sunlight or 
moonlight limitations. However, flux evaluations are possible imposing some 
restrictive hypoteses only: uniformity in space and time, and the mass index 
constant and equal to a mean value (Foschini, 1997). Radar data allow us 
to make a meteoroid streams map, different from that obtained with visual 
observations. Nevertheless, in both cases, it is possible to deduce that catas- 
trophic impact is a rare event, even if when there is a meteor storm (Beech 
and Brown, 1994; Beech et ai, 1995; Foschini and Cevolani, 1997). When 
micrometeoroids are considered the impact is not so rare: for example, there 
is 41% of impact probability for a space station (1000 m 2 area, 1 hour ex- 
posure), with a meteoroid with mass equal or greater than 10~ 8 kg, if there 
will be a storm like 1966 Leonids (Foschini and Cevolani, 1997). Such a flux 

3 For example, when the Halley comet reach the Sun, about the 3% of its surface is 
active. 

4 Some images of these structures can be found in: Beech and Brown (1994), Hughes 
(1995), Arter and Williams (1997). 
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can make serious damages to mechanical structures, particularly solar arrays 
and antennas, that can not be obviously shielded. 

4 Impact-generated plasmas 

The Olympus satellite experience, the post-flight analysis and the calcula- 
tions of impact probabilities impose a revision of potential space dangers. 
The mechanical impact do not seems to be a risk, as shown by several cases 
like the Hubble Space Telescope (Herbert and McDonnell, 1997) and the Mir 
space station (Christiansen et al., 1997). However, the Olympus failure is 
a paradigmatic example: in that case, the impact with a Perseid meteoroid 
could have generated electrical failures, leading to a chain reaction which cul- 
minated with an early end of the mission. According to Caswell et al. (1995), 
a gyro motor stopped, probably owing to a lack of power, and the satellite 
lost the reference. Following manoeuvres in order to acquire a new reference 
(Emergency Sun Acquisition) failed, probably owing to a short circuit in a 
capacitor of the emergency network. Even if there is not any certainty, it 
seems that the impact of a small meteoroid may have generated a plasma 
triggering discharge of charged surfaces, entering the grounded spacecraft via 
the umbilical. 

After the Olympus end-of-life anomaly, other authors looked to impact- 
produced plasma, rather than the impact itself (Brown et al, 1996; McDon- 
nell et al, 1997a). It is well known that, during an hypervelocity impact, a 
fraction of the projectile and target materials is evaporated and even ionized 
(Fechtig et al, 1978). A plasma cloud is then created almost instantaneously 
after the impact and expands into the surrounding vacuum. McDonnell et al. 
(1997a) find an empirical formula for evaluation of charge production during 
an hypervelocity impact. This equation, rearranged in order to emphasize 
projectiles dimensions and densities, can be written as: 



where 5 is the meteoroid density [kg/m 3 ], r is the meteoroid radius [m] and 
V its speed [km/s]. 

Because of the energy range, the plasma production is related to chem- 
ical composition of meteoroid. Cometary streams, richer of low ionization 




(5) 
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potential elements, will be more dangerous than other. The Leonid mete- 
oroid stream results to be the most dangerous stream, even during normal 
condition (McDonnell et al, 1997a). 

The impact-produced plasma can disturb the satellite in several ways: 
if directly injected into circuits can destroy part of the onboard electronics 
(McDonnell et al, 1997a); thermal forces can magnetize the neighbourhood 
of craters (Cerroni and Martelli, 1982); electromagnetic radiations emitted 
from the plasma can disturb several resources or scientific experiments on 
the satellite (Foschini, 1998). Even if satellites are actually submitted to 
several procedures for electromagnetic compatibility (EMC), meteoroid im- 
pacts call for new studies on these arguments. For example, the plasma- 
generated charge can deposit on near surfaces and, subsequently discharges 
to mass. The pulse shape will depend on electric characteristics (resistance, 
inductivity, capacity) of employed materials. Moreover, the pulse can dis- 
turb the onboard electronics, mainly in four ways (Audone, 1993; Foschini 
and Gallerani, 1993): 

1. the discharge can be directly injected into a circuit; 

2. the discharge can hit a nearly surface and disturb a circuit by a sec- 
ondary discharge; 

3. capacitive coupling between the discharge electric field and the circuit; 

4. inductive coupling between the discharge magnetic field and the circuit. 

If the first two modes are localized, and then depend on the impact 
place, the third and fourth modes can disturb distant components. How- 
ever, these coupling effects are strongly non-linear and depending on circuit 
layout. Thus, more detailed studies can be made on specific satellite only. 

5 Conclusions 

The threat from meteoroids must be revised, taking into account experi- 
ences such as Olympus and experimental studies (McDonnell et al, 1997a). 
Mechanical damages are localized, sporadically hit important parts and the 
catastrophic impact is an event still rare. On the other hand, if the plasma 
charge and current production are considered, then the risk increase and me- 
teoroid streams, particularly those composed with low ionization potential 
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elements, can be dangerous even during normal conditions. More studies 
about electromagnetic interferences from impact-produced plasmas are re- 
quired. 
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